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Abstract 

A class of meson decay modes sensitive to only one quark topology at leading 
Gp order (annihilation of valence quarks through a W) is described. No 
experimental observations, nor even upper limits, have been reported for these 
decays. This work presents a simpleminded (order-of-magnitude) calculation 
of their branching fractions, and compares to results of previous calculations 
where available. Although rare, one of these modes (D^ P^l) niight 
already be observable at charm experiments, two others (D"^ K*^'y, 
Dg'^'j) should appear at the B factories, and the rest at hadron colliders. 
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I. INTRODUCTION AND MOTIVATION 



As data on charm and beauty decays accumulates at CLEO, BES, and the Fermilab and 
LEP experiments, and soon at the B factories BABAR and BELLE, it becomes possible to 
study ever more rare decays in search of interesting and exotic physics. This work suggests 
examining a largely overlooked class, electromagnetic decays of a charged meson mediated 
by the weak annihilation of the meson. Such decays are indeed rare, with branching ratios 
(BRs) suppressed by a, not to mention the difficulty of forcing a camel (the whole meson 
wavefunction) through the eye of a needle (the pointlike W vertex multiplied by CKM 
factors), but their BRs need not be so tiny as one might think. 

Let us focus upon processes with only one meson in the final state. Such decays are 
especially interesting because they proceed through only one weak decay topology at leading 
order in Gp if i) the fiavors of all the valence quarks in both initial and final states are 
distinct, and ii) the initial and final quark have the same electric charge (and similarly for 
the antiquarks). This is the s-channel annihilation topology presented in Fig. |I]. Note that 
the first condition requires each valence quark to terminate on a fiavor- changing vertex, 
while the second eliminates the possibility of a t-channel W exchange. The photon may 
be attached to any charged line, although of course couplings to the lighter constituents 
are favored. At 0{G^p) corrections enter through diagrams with a penguin loop on each of 
the quark and antiquark lines (Fig. ^a), crossed-box diagrams (Fig. |^), and the diagram 
with the photon coupled to the W, none of which is expected to be very large. In addition, 
one may go beyond the valence diagrams and describe the weak process including its short- 
distance QCD corrections in terms of operators mixed through evolution of the corresponding 
anomalous dimension matrix,[| but we do not perform this refinement in this work. 

The interest in such decays lies partly in the simplicity of the weak topology and sensi- 
tivity to a number of hard-to-isolate CKM elements (as well as strong and electromagnetic 
matrix elements), and partly in the simplicity of the two-body final state. Indeed, the hard, 
monochromatic final-state photon should prove an exceptionally unambiguous experimental 
signal for these decays. It should also be pointed out that these modes would represent 
the first electromagnetic decays observed for the charged 0~ mesons D~^, D^, (except 
for the famous penguin mode B^ K*^--/), or B^. Table H presents decays representing 
the 6 possible weak annihilation fiavor assignments obeying our constraints, along with the 
CKM factors in the amplitude and their behavior in powers of the Wolfenstein parameter 
A 0.22. 

To our knowledge, the Cabibbo-unsuppressed mode p+7 was first studied in Ref. 

0] via quark model, then through pole and vector meson dominance (VMD) calculations 

light-cone and effective field theory techniques 0. Estimates for BR xlO^ vary from 
2.1 to 80 (see Ref. 0] for a summary). The present calculation, which for simplicity 
only takes into account one pole diagram in the language of Refs. gives 8 x 10~^. 

The double Cabibbo-suppressed mode — > K*^j was also studied in Refs. [QH, with 
BR results ranging from 3-30 x 10^^; we obtain 6 x 10^^. Encouraged by this consistency 



^Here we refer to mixing of the usual four-fermi operator with its Fierz reordering. One finds ||T],^ 
that the coefficient of the original operator can be enhanced by 20% or more. 
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with the more elaborate calculations, we apply our simple picture to the and modes. 
While D'^ — >■ p"'"7 is exciting because it might already appear in charm experiment data, 
— K*~^'y is interesting because it exhibits a neturinoless decay sensitive to |Ks|. 

The modes and D*+7 (collectively, D'^^^y), were suggested [|] as probes of 

\Vub\ and were estimated to have BRs of approximately 2 x 10^'' and 7 x 10"^, respectively. 
From Table |I| we see that the B~^ decays suffer the worst CKM suppressions. A number of 
the theoretical uncertainties associated with these estimates can be eliminated if the on-shell 
photon is replaced by an pair and the invariant mass of the virtual photon is 
used to define an operator product expansion. The price one pays for this improvement 
is an extra factor of a, so that such decays are estimated to have BRs of a few times 
IQ-w 10^^^, depending upon the level of Cabibbo suppression. The t-channel exchange 
processes mentioned above, which have neutral initial- and final-state mesons, are discussed 
in Ref. [0], and yield similar BRs. Most of these processes are too rare to be seen in 
appreciable numbers at the B factories (with combined yields of some millions of B^B~ 
pairs per year [|1T|), but may be observable at hadron collider experiments. 

The Be modes appear to have been studied only using the light-cone approach, in Ref. . 
We obtain BR results several orders of magnitude larger than theirs, and comment on this 
discrepancy below. 

The rest of this paper is organized as follows: In Sec. we present a very simple- 
minded calculation of the rate for these processes and a list of approximations used, while 
Sec. Ill presents numerical results, outlines experimental prospects for the observation of 
these modes, lists potential theoretical improvements, and concludes. 



II. CALCULATION 

The calculation presented in this section is very simplistic, in that it relies on a number of 
substantial approximations made explicit below. However, it is significant not in providing 
an exact determination of widths and BRs, but in obtaining the order of magnitude of these 
quantities as an estimate for experimenters searching for signals of these modes, and as a 
point of comparison for theorists performing subsequent, more refined calculations. 

The generic process we consider is M — >■ P ^ \^7, where M is the massive initial 0~ 
state, P is a lighter virtual 0~ meson with the flavor quantum numbers of the final state, 
and V is the final-state 1~ meson, as depicted in Fig. |^. For the charm decays, we have 
commented that consistency with the more elaborate calculations in Refs. indicates 
that this elementary ansatz seems to capture the essential order-of-magnitude physics. This 
approach avoids the danger of large cancellations between competing diagrams, but also 
runs the risk of missing important contributions in some cases. In modeling the decay this 
way, we make the following assumptions: 

1. Photon emission from M is neglected, so the process M — >■ M*7 — > V'^ is not included, 
as was done in Refs. Indeed, the D(^s)D^s)1 couplings have only measured upper 

bounds. In the charm case, Refs. used D(^s)D*g-^j input from previous theoretical 
calculations. In Ref. §, where M = B+, M* = B*+, P = P)+ , and V = D*+, the 
MM*^ and PV'j couplings were related through heavy quark spin-flavor symmetry 
(HQS) to that of DD*'j, and both diagrams were included. However, in the current 
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case with, for example, V = K*^ or p+, this is no longer an acceptable approximation, 
and we include only the photon coupling to the lighter mesons. This assumption likely 
leads to an underestimate of the BR, but probably not an exceptionally large one: 
The K*^ and electromagnetic widths are 50 ± 5 and 68 ± 7 keV, respectively, while 
that of the D*'^ is less than 4.2 keV;^ moreover, the calculations of Refs. show 
the D{^s)D*(^s)l amplitudes to be about 1/4 as large as those for PV'~^. Conversely, 
the light-cone calculations include only the photon coupling to M, arguing that 
couplings to V are suppressed by light quark masses. 

2. Complete factorization with vacuum insertion approximation is assumed for the weak 
vertex. The annihilation of M and the creation of P are assumed to occur at a 
single point. This approximation neglects both the short-distance QCD corrections 
as mentioned above, and long-distance hadronic contributions. As an example of the 
latter, the initial weak vertex may, rather than annihilating the initial meson, produce 
a quark (or antiquark) q that is the antiparticle of one of the meson valence quarks, 
and this four-quark intermediate state propagates for some time before qq annihilation 
occurs. Specifically, processes like [K'^tc^ or K^tt~^) — > K*^^, which require 
an understanding of final-state interactions, are not included, nor are the significant 
VMD diagrams such as K*^p^ K*'^'-/, where the p couples resonantly to a 
photon; nevertheless, the CKM coefficient of all these diagrams is the same. 

3. The intermediate state P is assumed to be the lightest pseudoscalar with the same 
flavor quantum numbers as the flnal-state V. Certainly many other resonant as well as 
multiparticle states with total angular momentum (that of M) can couple the weak 
vertex to Vj; however, the parity-violating couplings to O"*" states are neglected here. 
The present assumption is made partly because data exists on the PYj coupling from 
the observed decay V P7, and partly because the lightest pseudoscalar P among 
all possible intermediates presumably has the one of the largest couplings to V'j due 
to a relatively large wavefunction overlap. In any case, this approximation leads to an 
underestimate of the correct BR. 

4. In comparing the virtual process P — > V'-f to the on-shell V P7, one relates the 
single (magnetic) form factor C(g^) at a virtuality of = — mp to that at = 0. 
We take them numerically equal, although this tends to overestimate the rate, since 
form factors tend to fall off away from = 0. Nevertheless, we indicate this ratio 
explicitly in the flnal expression, see Eq. (|2.4|) below. 



Given these assumptions, the calculation of the rate is a simple matter. The weak mixing 
vertex, mediated by an operator Ow, in vacuum insertion approximation is given by 



{PipM)\Ow\MipM)) 

V2 



-t^VpVMB {p{pm) Qp7^(l-75)gp 0)(0 g^,7M(l-75)QM M{pm] 



^Nonetheless, one should note that a small propagator denominator {niy — m\.jt) or different a 
Ught quark charge in M* can enhance the importance of such couplings in the full width. 
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V2 

-t^VpVMfpfuM'B, 



mPm I 



(2.1) 



where ttim is now abbreviated as M, the valence structure of P is QpQp, that of M is Q^^q^-j, 
and Vm,p are the CKM parameters associated with the annihilation of M and creation of 
P, respectively. We also allow for a coefficient B parameterizing the incompleteness of the 
vacuum saturation approximation, as in BB mixing, but set it to unity in our numerical 
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estimates. The intermediate state P provides the simple propagator i/{p' 
photon vertex is extracted from the decay Ty^p-y, which has invariant amplitude 



m. 



M = C{p 



m. 



The 



(2.2) 



where the Lorentz coupling is determined in part by the 0~ quantum number of P. C 
is a magnetic form factor, and simply becomes a transition magnetic moment when its 
momentum transfer argument pp — mp is set to zero in the on-shell case. The rate obtained 
from this amplitude is 
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(2.3) 



The full rate for M — > Vj also uses C^, but now has the argument — mp. For our 
numerical estimates, we assume that C does not change dramatically over this range, and 
use data on Ty^p^y to eliminate C from the expression for TM^v-y] however, since this is 
certainly a contentious approximation, we formally retain the ratio of C at the two different 
argument values in the full expression for the width. Putting this together, one obtains our 
central result: 



r(M ^ 1/7) 



-Gp I VmVpI^ fljfpB'^ Ty- 



'C{M^ - ml) 




(2.4) 



Here, the cubed mass factors are nothing more than the ratio of E'^ for M V'y to that 
for V P-f. 

One may compare Eq. ( p.4|) to Eq. (17) of Ref. |^ for B^ ^{s)' which in the current 
notation reads 



27 



r(M ^ V^) = \VmVp\^ ft,B^ Tv'^p 



'C{M^ - ml) 



X 



mvM{M - mv){M + myf 
{my, — mp,)^ 



m 

. M 



(2.5) 



where the 0~, 1^ pair P', V are related to P, V by HQS: The primed mesons are introduced 
when data for Fy^p^ is unavailable. In Ref. P', V were , D*^ and the experimental 
upper bounds for D*^ D^'j were used; the current calculation uses only information 
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from the unprimed mesons directly, so P',V' P,V here. All of the differences between 
Eqs. ( p.4|) and (|2.5| ) can be accounted for by the assumptions of HQS: First, the magnetic 
moment form factors for all heavy mesons were assumed the same, except for a trivial 
coefficient due to the electric charge Qq(m) of the lighter quark q in the meson m to which 
the photon couples. Since both the MM*7 and PV'j couplings were included in the HQS 
calculation, an extra enhancement of {Qq{M) + QqiV)Y /Ql(v)y ^ factor of 9, appeared in 
Ref 0. Next, in HQS one has fljM = f^mp, and -mp = 0(1/M), M being the heavy 
quark mass. The remaining differences arise from the fact that fields containing heavy 
quarks in HQS are nearly static, even if the heavy quark changes fiavor. This leads one 
to adopt the normalization of HQS states of 1 rather than 2M particles per unit volume, 
as well as introduce propagators linear rather than quadratic in particle masses, and these 
modifications often lead to effective lowest-order substitutions (in the current notation) such 
as (M + mv)/2 M. Equations ( ^.41 ) and ( |2.5| ) are related by the application of these 
properties. 



III. RESULTS AND CONCLUSIONS 

We employ Eq. (|2.4|) to obtain BR estimates for the 6 modes exhibited in Table H We 
use standard Review of Particle Physics (RPP) |]I2[ values for masses, decay constants, 
CKM elements, and lifetimes whenever possible, with the following exceptions: We take 
\Vub\ = 3 X 10-3, = 170 MeV, and fo = fo, = 200 MeV. Following recent experiments. 



we use the E791 value p| td, = 0.518 ± 0.014 ± 0.007 psec, and the CDF values 114 
= 6.40 ± 0.39 ± 0.13 GeV, = OA6to\l ± 0.03 psec. OPAL |15| and ALEPH p 



have also reported a few Be candidate events, with mass values consistent with Ref. []I4 
We also use the most recent lattice determinations of which combined give 425 ±11 
MeV. As mentioned previously, we take B"^ = 1 and C^(M^ — mp)/C^(0) = 1. 

The resulting BRs are exhibited in Table ||, along with the energies of the final-state 
photon. We see that the B^ modes give BRs in agreement with Ref. despite a very 
different calculation, and similarly for the decays in comparison with Refs. as 
promised in Sec. |. A few — > K*^'-) and B^ —>■ D*~^'-f should appear each year at the B 
factories, with the exact number depending upon the correct value of the mantissa in our 
estimate. The mode D'^ — > p'^j might even be observable right now at charm experiments 
such as at BES (or possibly CLEO) if our estimate is low by a factor of a few, or the upper 
bounds of the estimates in Refs. |^,^ are correct, based on limits in the RPP [|T2[; in any 
case, existing experiments can place a meaningful upper limit on its BR. 

However, in contrast to Ref. we find that the B^ modes are rare but not exceptionally 
so; they find BR{Bc p^-f) = 8.3 x 10"^ and BR(5c K*+-f) = 5.3 x 10"^ This 0(10^) 
discrepancy might be explainable if C(0)/C(mp — m^) ^ 30, but then the order-of-magnitude 
agreement for D(^s) and B decays becomes a mystery. Furthermore, even if the form factor 
C falls off this dramatically, other longer-distance mechanisms (see point 2 in Sec. H) would 
likely step in to maintain the rate. Of course, since the B factories are not designed to 
produce B^s, the possible observation of these modes must necessarily wait for the upcoming 
hadron collider experiments at the LHC or Tevatron. 

Although much physics is neglected in this simple calculation, our estimates show that 
weak annihilation decays may be observed in the near future. They are attractive from both 
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the experimental and theoretical perspective. Improvement of the theoretical calculation 
essentially amounts to improving on the four assumptions made in Sec. |l[ Lifting the 
first requires new data for the heavy coupling in M* M7, particulary positive data for 
D*^ D^'^. Improving the vacuum insertion approximation can be accomplished as in BB 
mixing, with lattice or model calculations, and including the short- distance QCD corrections 
is a straightforward matter. Many of the neglected long-distance corrections such as VMD 
diagrams have been considered in Refs. but one must take care with their relative 

phases, while final-state interactions must still be taken into account. As for the remaining 
two assumptions, one can be freed of both the lowest-resonance dominance and constant form 
factor assumptions by either including other intermediate channels explicitly, or carrying out 
light-cone or inclusive quark model calculations. 

From the experimental side, it would be interesting to see what direct bounds can be 
placed on these modes at the current time, in anticipation of their eventual observation. Once 
observed, the weak annihilation modes will present an interesting probe of CKM elements, 
electromagnetic transitions, and meson wavefunctions. 
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TABLES 



Valence structure 


Decay mode 


CKM Elements 


bu — > CS7 


B+ - 


. Z)*+7 




bu cdj 


B+ - 


^ D*+7 




be — > us'y 


Bt- 


> i^*+7 




be — > duj 




p+7 


nVud - A^ 


cd — > its7 


D+ - 


> iv:*+7 


KdKs ~ A2 


cs — > udj 


Dt- 




KsKd ~ AO 



TABLE L Flavor structure and mesonic decay modes of weak annihilation electromagnetic 
decays. The CKM coefficient for each process is accompanied by its magnitude in powers of 



Wolfenstein A « 0.2. 



Decay mode 


BR (est.) 


Photon Energy (GeV) 


B+ ^ Z):+7 


1 X 10"'^ 


2.22 


5+ ^ L)*+7 


7 X 10-9 


2.26 


B+ K*^-/ 


3 X 10^6 


3.14 


B+ ^ p+7 


3 X 10"^ 


3.15 


D+ A*+7 


6 X lO"'^ 


0.72 




8 X 10"^ 


0.83 



TABLE IL Estimates of branching ratios for weak annihilation decays using Eq. (2.4). Also 
included are energies of the monochromatic photon. 
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FIGURES 



FIG. 1. The s-channel weak annihilation topology. 



(b) 





FIG. 2. 0{G^) corrections to Fig. ||: (a) di-penguin diagram; (5) crossed-box diagram. 



M 




FIG. 3. Meson decay diagram for the process M 
vor-changing vertex of Fig. |l| and its corrections. 



P V^. The blob represents the fla- 
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